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Abstract. Loops are subgraphs responsible for the multiplicity of paths going from 
one to another generic node in a given network. In this paper we present an analytic 
approach for the evaluation of the average number of loops in random scale-free 
networks valid at fixed number of nodes N and for any length L of the loops. We 
bring evidence that the most frequent loop size in a scale-free network of N nodes is of 
the order of N like in random regular graphs while small loops are more frequent when 
the second moment of the degree distribution diverges. In particular, we find that 
finite loops of sizes larger than a critical one almost surely pass from any node, thus 
casting some doubts on the validity of the random tree approximation for the solution 
of lattice models on these graphs. Moreover we show that Hamiltonian cycles are rare 
in random scale-free networks and may fail to appear if the power-law exponent of the 
degree distribution is close to 2 even for minimal connectivity kmin > 3. 
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The scale-free network structure has been found in a number of social, technological 
and biological networks as the skeleton of their interaction P3 121 E] ■ The main property 
of scale-free networks is to have a power-law degree distribution P{k) ~ k~"' and second 
diverging moment, i.e. 7 G (2,3]. To distinguish between different scale-free networks, 
recently, much attention has been devoted to network motifs jHEllE], i-e. subgraphs 
that recur with higher frequency than in maximally random graphs with the same degree 
distribution. Among those, the most simple types of subgraphs are loops [HIHlin], i-e. 
closed paths of various length that visit each node only once. Loops (or cycles) are 
interesting because they account for the multiplicity of paths between any two nodes. 
Therefore, they encode the redundant information in the network structure. Another 
discriminant aspect of real scale-free networks is the presence of degree correlations 
between linked nodes. Characteristic motifs in a graph and degree correlations are in 
many real graphs not independent phenomena but they depend on each other as it has 
been shown for small (up to maximal connectivity) size subgraphs in ^1 [TH] . 

Last but not least, it has been observed that the distribution of loop sizes is intimately 
connected with the thermal properties of lattice models defined on that graph [THl I17j . 
On the other hand, the analytic approach to these models relies on the assumption that 
locally a random graph can be considered to have a tree like structure [1211201121], i-e. 
that loops of finite size are rare. 

In this paper we present an analytic derivation of the average number of loops of 
any size in a random scale-free network. Our motivation is that the results on random 
networks provide a reference picture which often captures key intuition which extends 
to correlated networks (see e.g. ^lE])- In addition, it provides valuable information 
for the statistical mechanics of lattice model on random graphs |18| I17j. 

Let us first recall the classic results for regular random graphs, i.e. random graphs 
with fixed connectivity of the nodes fcj = c for each node i. A regular random graph 
contains a finite number of small loops of size L ^ log(A^), with average expected 
number 



and Poisson fluctuations around the mean. On the contrary for large loop sizes 
L ~ 0{N) the number of loops goes as 



where i = L/N and cr(£) is a function having the maximum at imax = c/{c + 1) whose 
expression can be found in the literature [inilllj. Regarding Hamilton cycles, i.e. loops 
that span the entire network L = N their expected number for a large regular random 
graph is diverging with the system size as long as c > 3. For c = 2 the average number 
of Hamilton cycles goes to zero as the system size diverges ^Hl • Coming to the scale-free 
network literature, Ref.0 analyzes the number of loops of any size on a pseudo-fractal 
scale-free graph and report the scaling behavior 
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with L* = A^V{7-i). ^gg^i^ 

were presented on Hamilton cycles, to our knowledge, so 

far. 

In the following we characterize the statistics of loops in random scale-free networks. 
We find a larger number of small loops with respect to regular random graphs. In 
particular, we compute the expected number of loops of a given size passing through a 
node and find that when 7 G (2, 3) this number diverges with the network size, beyond 
a finite loop size. This raises some doubts on the solution of lattice models on these 
graphs based on the local tree approximation [ini EDI EII • We also find that loops have a 
characteristic size L* ~ A^. In other words, our results are consistent with the scaling Q 
with L* ^ N and with Eq. Q for regular graphs. This suggests that the result of Ref. 
[Hj crucially depends on the peculiar correlations of the ensemble they consider. Special 
attention will be given to Hamilton cycles that in networks with a small 7 exponent can 
fail to exist unless the lower cutoff of the distribution is large enough. 

There are different ensembles of random networks one can consider. The classical 
one follows the prescription of Molloy and Reedj21I: First, to each node i of the network 
is assigned a connectivity fc, drawn from the chosen probability distribution, secondly 
edges are randomly matched. This ensemble indeed generates networks of given degree 
distribution but it may yields networks with multiply occupied links. More precisely, the 
distribution of the links between two nodes of connectivity fcj, kj is a Poisson variable 
with mean kikj/{cN), where henceforth c = (k) will denote the average connectivity. 
Hence the probability of no multiply occupied links is 



where the right hand side refers to a scale free random graph with degree distribution 
p(k) = ak-'^ with k G [m, K]. Taking a structural cutoff K ~ iV^/^ [23 EEI, we 
conclude that double links appear with probability one for 7 G (2,3] as iV — > cx) in the 
Molloy-Reed ensemble. When counting loops of a network this effect becomes relevant 
and undesirable. Thus we consider also another ensemble where double links cannot 
appear: the static fitness network I23j. In the fitness ensemble nodes are assigned 
a random variable (fitness) q drawn from a p{q) distribution function and every couple 
of nodes is linked with a probability depending on the fitness of the considered nodes 
p{q, q'). When p{q) is power-law distributed and p(g, q') = -^j;^ the resulting graph is a 
random scale-free graph characterized by the same exponent of the fitness distribution. 
In these graph the connectivity of every node is a Poisson variable with expected value 
{k{q)) = q. This ensemble doesn't allow for networks with double links but instead may 
give rise to networks with isolated nodes {ki = 0) or to nodes connected with a single 
link {ki = 1) to the others. The presence of such nodes rules out the possibility to find 
Hamilton cycles, hence we shall take this effect into account when discussing Hamilton 
cycles. 

Consequently, the Molloy-Reed ensemble and the fitness ensemble are not equivalent 
and have intrinsic properties that could perturb sensitively the counting of the number 
of loops. In order to understand the dependence on the details of how graphs are 
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generated in the following we are going to study the expected number of loops in the 
two ensembles. 

1. Loops in the fitness ensemble 

The prescription of Ref.j22| to generate a class of random scale-free networks with 
exponent 7 is the following: i) assign to each node i of the graph a hidden continuous 
variable distributed according to a scale-free distribution p(g) = poq^'^ for q G [m, Q] 
with po = (7 — 1)/ {ni^^'^ — Q^^^) the normalizing constant. Then ii) each pair of nodes 
with hidden variables g, g' are linked with probability qq'/{cN), where c = (g) is the 
expected value of g. For 7 < 3 the cutoff Q ~ N^^'^ is needed to keep the linking 

probability smaller than one, i.e. Q'^/{cN) < 1 while for 7 > 3 the cutoff is the natural 

1 

one Q ~ iVT-i . By construction the expected value of the connectivity of a node with 
hidden variable g is (/c|g) = g and there are no multiple connections between nodes. 
Notice that the average connectivity of the graph is given by 

7 - 1 

(k) = {q) = c^ -m (5) 

7-2 

in the limit N 00. 

A loop of size L is an ordered set of distinct nodes {Ij, . . . , z^}. For each choice of 

the nodes, the probability that they are connected in a loop is 

2 

qiiqi2 qi2^i'i qiL^h tt 'iii 

Nc Nc ivT ~ Nc' 

The total number of possible loops joining these L nodes, in any possible way is L\/ (2L) 
where the factor 2L comes from the fact that the initial node of the loop can be chosen 
in L ways and that there are two orientations. In order to count loops, let us lump 
together nodes with hidden variable g, G [g, g + Ag), where Ag is a small interval of g. 
In each interval of g there are Ng ~ NP{q)Aq nodes of the network. For each choice 
of the L nodes, let Ug be the number of nodes with g^^ G [g, g + Ag). Then the average 
number of loops of size L in the graph is given by the number of ways we can choose {ug} 
nodes multiplied by the probability that these nodes are connected in all distinguishable 
orderings. Consequently we have 

where the sum is extended over all {rig} such that J2q = L. Introducing this constraint 
with a delta function and using its integral representation, we find 

J^^^ ^ f°° ^^iLx+N {log[lWe-'V (Nc)]) ^ 

2L J-00 271 

Notice that in Eq. ((7j) one can safely take the limit Ag —>■ and that the average over 
the -P(g) distribution is taken assuming that we focus on the limit — 00. In what 
follows, we will evaluate Eq. ((Tj) in different ranges of L < A^ in the limit A^ ^ 00. 
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1.1. Small loops 

For L finite but large, the integral in Eq. ((Zj) is dominated by values x ~ —iz* where 

. M f 1 _ M£ + ] (8) 

- Lc [ {qY N I ^ ' 

where we have neglected all terms beyond the first leading correction when N oo. 
The argument of the exponential in Eq. ((Tj) can be expanded around x ~ —iz* yielding 
(log [1 + q'^e-'^'/cN]) + Lix ~ z* -^x- iz*Y + 0{x - iz*f'\ . Hence the 

integral can be estimated by saddle point for L large. Using the asymptotic expression 
L\ ~ y/27cLL^e~^ , we find to leading order 

L 

(9) 




This approximation is valid as long as the leading correction in Eq. (jHl) is small. Using 

that (g") = Po(Q""^"^^ - - 7 + 1) for 7 ^ n and that Q ~ mm{N^/^, N^) 

we find that the expression above for A/l holds when 

N 7 > 5 



L < ~ <! 7 > 3 (10) 

^ [ iV(3-7)/2 2 < 7 < 3 

with logarithmic corrections for 7 = 3 and 5. Note that strictly speaking the expansion 
(jH)) is converging only for N{q^)/L > A^, i.e. L < A^(3-7)/2 for 2 < 7 < 3 and 
L <ti ]\[i'y-'^)/i'y-'^) for 7 > 3. Nevertheless Eq.Q remains valid in the limits (fTIH) as an 
asymptotic expansion. For 7 > 3 we obtain a result very similar to Eq. (P) for regular 
graphs. On the contrary, for 2 < 7 < 3 we have (g^) ~ aN^^~"'^^'^, with a a constant, 
hence the number of finite loops 

Ml^^I-Yn"^^ (11) 



2L\c 
diverges as 00. 

1.2. Intermediate loop sizes and the most frequent loops 
It is convenient, at this point, to write Eq. ((Tj) as 

^ r _^^NfiNce^^,L/N) (^2) 

J-oD 2 V ttL 

where we have used Stirling's approximation and 

f{y, i) = (log [l + q^y] ) + i log{iy/c) - L (13) 

The integral can be computed by saddle point method, deforming the contour 
of integration so as to pass from the point where / is stationary. The condition 
9y/(l/,^) = yields 
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Figure 1. The behavior of £max as a function of 7 for m = 1,2,3 and N = 10^ 
for fitness ensemble (solution of the Eq. H17|l in solid lines) and for the Molloy-Reed 
ensemble (solution of the Eq. H32|l in dashed lines). 



Let y*{i) be the value of y which solves this equation. We can expand f{Nce^ 
the corresponding (complex) value x* of x 



/(iVce", e) = f{y\ i) - ^—^ix - x*f + . . . 



2 dy^ 

V* ' 



2 \ [g^ + y* 



ix - x* f + 



As long as 



> 1 



we can neglect higher order terms. This yields the leading behavior 



Mr 



n-1/2 



LNy* 



[q^ + y* 



around 



(15) 



(16) 



The number of loops A/^ takes its maximum for loops of length L = Ni^^^ where 



The solution £max is plotted in Fig. [T] against 7 for scale-free graphs and different values 
of m, for = 10^. Notice that as 7 — 2"*", the size of most probable loops vanishes as 
^max ~ 7 — 2. Around the maximum, A/l takes a form similar to that for regular graphs 
(see Eq. 12)) which is consistent with the scaling form Eq. with L* ~ A^. 

For scale-free random graphs with 2 < 7 < 3 there is an intermediate range of loop 
sizes L ~ A^(3-7)/2 -y^^j^ich is related to the solutions with y* = ijlN with yU > 1. More 
precisely, we find that for loops of size L = xilA-^^'^"^'^^'^ have 



Ml 



G(/i) 



,L[log(ML/c)-l+//(A.)] 
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with G{fi) a function of /i and 



x(/x) = (7-iw-^E(-ir: 



ra=0 



3 - 7 + 2n 



^1 - 7 + 2n) 

Notice that A/l does not satisfy a simple scahng form such as Eq. (jH]) in this intermediate 
region. 

1.3. Hamilton cycles 

From Eq. (jH)) we can easily calculate the number of Hamilton cycles. Indeed for L = N 
we have the asymptotic behavior 



This is the expected number of Hamiltonian cycles over all the networks of the fitness 
ensemble including networks with nodes of low degree fcj = 0, 1, which by definition 
cannot have an Hamilton cycle. It seems a sensible thing to compute the number of 
Hamilton cycles in networks with a minimal degree connectivity grater than 3, i.e. 
ki >= 3. In fact it is well known that for a regular random graph of connectivity c = 2 
the expected number of Hamilton cycles goes to zero in the N ^ oo limit whereas 
regular graphs are Hamiltonian when c > 3. Taking this as a reference result, we 
normalize A/at by the probability vr that all the nodes have at least 3 connections. Since 
the connectivity of each node in the fitness network is aPoisson variable with expected 
value q the probability that all the nodes have connectivity /c > 3 is simply given by 
vr = e^^(™'^), where 

A(m, 7) = (log(l - (1 + g + q'me-")) (20) 

In the limit N (yo we find 

_log__log(^-— ^mj +-— ^-A(m,7). (21) 

This implies that if random scale- free graphs have Hamilton cycles only for m > 771^(7) 
where mdj) is the value of m for which Eq. (j2T| vanishes. Conversely, for m < rriclj) 
a random scale-free graphs has almost surely no Hamilton cycle. 

In Fig. 121 we report the critical value mcij) as a function of 7. Notice that 
rric ~ 1/(7 — 2) — s> 00 as 7 — s> 2^. Consequently, if we consider only the networks 
of the nesemble with kmin > 3, we find that random graphs with 7 < 7^^^. = 2.16 .. . 
(where rnd'jpit) — 3) do not have Hamiltonian cycles. Considering that regular random 
graphs with k = c> 3 are Hamiltonian, this may seem a surprising result, at first sight. 
The basic intuition to explain this apparent paradox is that most paths pass through 
well connected nodes. Hence even if /cj > 3 it is very unlikely to have a path spanning 
the entire network which is not passing through the most connected nodes more than 
once. 
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1.4- Loops passing through a node 

In order to count of the number of loops of size L passing through a given node, with 
fitness value qi, we can repeat the previous outlined above, without taking the average 

3—7 

over Qi. For 7 < 3 and short loops sizes L -C this gives the expected number 

A/ite)-^T^f-V^ArV(^-'). (22) 



cN2L \cj 

Focusing on nodes with ~ A^", we find that the number of loops of size 

L>L^ = l + -^{l-2a) (23) 

3-7 

diverges with the network's size A^. For example, nodes with a finite qi have an infinite 
number of L = 5 loops passing through it in networks with 7 < 2.5 but at most a 
finite number of loops of size L = 3 if 7 > 2. The most connected nodes (a = 1/2) 
instead have an infinite number of loops of any size L > 3 passing through it. Notice 
that Lq 00 as 7 — > 3 in order to match the behavior Lq ~ logA^ of regular graphs. 
Conversely, in a finite graph of N nodes, only the large fitness nodes with 

g. > AT^- V(^-i) (24) 
belong to a significant number of loops of size L. 



2. Loops in the Molloy-Reed ensemble 

The counting of the number of loops in the MoUov-Reed (24j follows a procedure much 
similar to the one considered for the fitness ensemble nevertheless giving different 
results. To construct a Molloy-Reed network one proceed as follows: i) a, degree is 
assigned to each node of the network following the desired degree distribution with 
cutoff K ~ min(ArV2^Ar— ). De gree distributions which do not satisfy the parity of 
cN = J2i ki are disregarded; ii) the edges coming out of the nodes are randomly matched 
until all edges are connected. When this procedure ends with nodes having links to 
themselves (tadpoles), the whole network is rejected and the procedure is started anew. 

To calculate Ml in this ensemble first one has to count in how many ways it is 
possible to have a loop of size L in the network and weight the results with the fraction 
of possible networks in the ensemble which contains the loop. Let us first state that 
the total number of graphs in the Molloy-Reed ensemble is given by {cN — 1)!!. Indeed 
when constructing the network by linking cN unconnected edges one start by taking 
one edge at random and connecting it to one of the {cN — 1) possible connections. 
Then one proceed taking another edge and linking it to one of the remaining {cN — 3) 
possible connections thus giving rise of one of the {cN — 1)!! possible networks. By 
similar arguments one shows that the total number of networks containing a given loop 
of size L are {cN — 2L — 1)!!. On the other side the total number of loops of size L in 
the Molloy-Reed ensemble are given by the number of ways one can choose an ordered 
set of L nodes {Ij, . . . , ii} of connectivity {ki, /c2, . . . , /cl} and connect them on a loop. 
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As for the fitness network the total number of possible loops joining these L nodes, in 
any possible way is L\/{2L). The number of ways one can choose the edges coming out 
of the nodes to form the loop is given by 

T^LMki - 1). 

Consequently the average number of loops in the MoUoy-Reed ensemble will be given 
by 

= ^ E n , u'" , . (Kk - i)r w.,, (25) 

where Nk = NP{k) (uk) is the number of nodes with connectivity k present in 
the network (loop), K is the cutoff of the degree distribution and the sum over 
{uk} is restricted to {n^} such that J2k^k = L. Moreover we use the definition 
Wn,l = {cN — 2L — l)\\/{cN — If we use the Stirhng approximation for Wn,l 
we get the expression 

Wn,l ~ (cAr)-^e^^(') (26) 

with i = L/N and 

g(i) = ^(c-2i) log (^^^^ +L (27) 

Thus we get 

M.^Iiyu W { k(k-i) \ „„„ 

which except for the substitution k{k — 1) and the factor exp(iV5f(£)) is equivalent 
to the expression ((7j) of the average number of loops of size L in the fitness ensemble. 
Following the same steps as in the fitness ensemble, we get 

J\J-^ ^ 1"^ ^^iLx+N{log[l+k(k-l)e-^^/(Nc)])+Ng{e) ^29) 

2L J-oD 271 
with g{i) given by Eq. (j^Tj). 

2.1. Small loop size 

The number of small loops in the Molloy-Reed ensemble is given by 

where this approximation is valid asymptotically for loops sizes satisfying Eq.()10j). Note 
that as in the fitness ensemble for 7 G (2, 3) short loops diverge as A//, ~ . 
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2.2. Intermediate loops sizes 

For intermediate loops in the MoUoy-Reed ensemble a similar expression to Eq. ((1 
holds with 

fiv, i) = (log[l + kik - l)/y]) +nogiiy/c) - i + g{i). (31) 

The calculations of the average number of loops is very similar for the MoUoy-Reed and 
fitness ensemble, with a difference for the equation of the loops of maximal size which 
in the MoUoy-Reed ensemble satisfy 

/ k{k-l) ^^2^ 

\ C '^^max ~\~ k(^k l)^max / 

In Fig. ^we report the value of imax in the MoUoy-Reed networks as a function of 7 for 
different value of the minimal connectivity m for = 10^. 



2.3. Hamiltonian cycles 

Starting with expression (j^^ one can easily evaluate the expected number of 
Hamiltonian cycles in MoUoy-Reed networks. Indeed ioi L = N and c > 2 one can 
use the Stirling approximation to find the asymptotic behavior (A^ 00) 

1 log(Ar^) = {\og{k{k - l)/c)) + \{c-2) log(l - 2/c). (33) 

If we approximate k{k — 1) with k"^ which is possible close to 7 ^ 2 in the limit c — > 00 
we recover the same behavior as in the fitness ensemble: if the minimal connectivities m 
is smaller than the value mc(7) for which Eq. (jH^ vanishes, then a scale-free network 
is typically not Hamiltonian. In Fig. |21we report mc(7) for 2 < 7 < 3 and we confirm 
the behavior ~ 1/(7 — 2) for 7 — > 2. For example, we find that MoUoy-Reed random 
scale-free graphs with minimal connectivity m = 3 are typically not Hamiltonian if 
7 < 7j^j^ = 2.27. . .. As for the fitness ensemble, the intuition is that it is not possible 
to extract from a random-scale-free graphs a subgraph which is a regular random graph 
with fixed connectivity c > 3 if m < mc(7). 



2.4. Loops passing through a node 

To count the loops of size L passing through a given node of connectivity k we must fix 
it and choose other L — 1 nodes to form the loop. For short loops sizes and exponent 
7 < 3 this gives the expected number 

J^L{k) ^ ^^^^Ar,_, (34) 

with Ml ~ N^^. The same results derived for the fitness ensemble hold: There is a 
critical finite loop size Lo{k) such that there are infinitely many loops of size L > Lq 
passing through a given node of connectivity k. On the contrary, in a finite but large 
network of A^ nodes, loops of size L becomes significant for nodes of connectivity 

A:»Ari-^(^-i). ^35^ 
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Figure 2. Dependence of mc(7) on 7 in the fitness (solid line) and MoUoy-Reed 
ensemble (dashed line) in the limit N 00. Observe that for 7 < 7* we have md^) > 3 
in both ensembles. 



3. Numerical results 

We compare the analytic results derived so far with the direct count of the number 
of loops in a sample of computer generated random graphs in both the fitness and 
the MoUoy-Reed ensemble. This is important because Ml is a fluctuating quantity 
which takes exponentially large values. In other words, the analytic calculation of the 
expected number of loops may be dominated by (exponentially in A^) rare realization of 
graphs with an exponentially large number of loops. In this case the number of loops 
of a typical realization of a graph would differ from our estimate. We have chosen the 
fastest known algorithm for calculating the total number of loops exactly j2ZI as in Ref. 
|17j . This algorithm has a upper time bound of 0{NC) where C is the total number 
of loops in the network. The simulations performed in this way enable one only to 
consider small networks sizes < 50 and small m < 3 as the total number of loops 
in such graphs increases exponentially with the system size. Note that for such small 
sizes the degree distribution contains nodes of very similar degree since the upper cutoff 
is K ~ 6 for 7 = 2.1. Moreover in order to compare the direct counting with the 
analytical calculation, we have chosen a fixed degree (fitness) distribution A^^ = NP{k) 
to reduce fluctuations that become relevant for such small sizes. In Fig. IHlwe report the 
analytic prediction of the average number of loops of a given size in a fitness network of 
= 30 nodes. This results are compared with direct counting of the loops in computer 
realizations of these networks were data are averaged over 50 realizations. We found 
strong sample to sample fluctuations which we believe are responsible for the deviation 
from the analytical results. 

On the contrary, for the Molloy-Reed networks of same system size the direct count 
of loops is very close to the analytic prediction. Fig. HI reports the direct count of 
loops for Molloy-Reed networks [28] with A^ = 30 and several degree distributions and it 
compares it with the corresponding analytic prediction. 
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Figure 3. Number of loops for fitness networks of = 30 nodes and given distribution 
p{q) with TO = 3. The average is taken over 50 realizations. 
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Figure 4. Number of loops for MR networks without double links of N = 30 nodes and 
different 7. The direct count averaged on 10 realizations is compared to the analytic 
prediction for the same degree distribution (full lines). 

4. Conclusions 

In conclusion we have computed analytically the expected number of loops of any size 
in a scale-free network. We found that scale-free graphs have a very large number 
of small loops compared to regular random graphs. On the contrary we have shown 
that, also with a minimal connectivity kmin > 3 the expected number of Hamilton 
cycles can be zero in the N ^ 00 limit provided that 7 is sufficiently close to 2. The 
reason for this is that paths connecting many nodes need to pass frequently on nodes 
with high connectivity. Put differently, it is not possible to embed a regular graph of 
connectivity c > 3, which would have an Hamilton cycle, in scale free networks if 7 
is too small, even if all nodes have ki > c. In the intermediate region of relatively 
large loops we found that the expected number of loops attains its maximum for loops 
of size L N. These results are derived both in the fitness and in the MoUoy-Reed 
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ensembles. While the generic picture is the same, the results in the two ensembles differ 
quantitatively highlighting that the loop size distribution is somewhat sensitive to the 
precise prescription for drawing random graphs. Moreover we have checked the results 
with direct counting of computer generated scale-free networks belonging to the two 
ensembles. It would be desirable to derive similar results for ensembles of correlated 
scale-free networks. 
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